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Abstract 

It is shown that the partition function of the 2d Ising model on the dual finite lattice 
with periodical boundary conditions is expressed through some specific combination of 
the partition functions of the model on the torus with corresponding boundary condi- 
tions. The generalization of the duality relations for the nonhomogeneous case is given. 
These relations are proved for the weakly nonhomogeneous distribution of the coupling 
constants for the finite lattice of arbitrary sizes. Using the duality relations for the non- 
homogeneous Ising model, we obtain the duality relations for the two-point correlation 
function on the torus, the 2d Ising model with magnetic fields applied to the boundaries 
and the 2d Ising model with free, fixed and mixed boundary conditions. 
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1. Introduction 



The duality relation for the two-dimensional Ising model was discovered by Kramers and 
Wannier [1] in 1941 year. They established the correspondence between the partition function 
of the model in low-temperature phase and the partition function in high-temperature phase 

(smh2K))~ N / 2 Z(K) = (smh2K)- N / 2 Z(K) (1) 

sinh2iT • sinh2X = 1. 

Using this self-duality property, in [1] the critical temperature in the 2d Ising model was 
determined before the Onsager's exact solution [2]. 

In paper [3] Kramers- Wannier duality relation (1) was generalized to the nonhomogeneous 
case (the coupling constants are arbitrary functions of lattice site coordinates) 

n(sinh2^(f))- 1/4 Z[^(f)] = n(smh2^(r))- 1/4 Z[^(r)], (2) 

r,i r >* 

sinh2iri(r) • sinh2^ 2 (r) = 1, 

sinh 2X2(7-) • sinh2^i(f) = 1. (3) 

Here r, f and Ki(r), Ki(r) are coordinates and coupling constants on the initial and dual 
lattices respectively (we will define them in the following section). Since Kadanoff-Ceva 
anzats (2) defines connection between functionals but no functions, this relation is very use- 
full for analysis of the thermodymamic phases of the model. So, for example, this anzats 
allows to correctly define disorder parameter /x, to obtain the duality relation connecting 
correlation functions on the initial and dual lattices, to define "mixed" correlation functions 
(cr(rj) . . . <7{rj),n{r k ) . . . fi(ri)) and e.s. 

As was already mentioned in [1,3] relations (1) and (2) can not be understood as literal. So, 
for example, using method of comparing of high- and low-temperature expansions for deriving 
of duality relation (1), it is hard to take into account and to compare the graphs, which include 
spins on the boundaries (in particularly, graphs which contain cycles wrapping up the torus). 
In fact (1) is correct in the thermodynamic limit (for the specific free energy). However for 
the nonhomogeneous case the procedure of thermodynamic limit is rather ambiguous. In any 
case it is usefull to have exact relations (in contrast to (1) and (2)) connecting the partition 
functions on the initial and dual finite lattices. This is aim of our paper. 

In the section 2 we introduce denotions and define the representation of the partition 
function in the form of the functional Grassmann integral. In section 3 we derive exact 
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duality relation for the homogeneous Ising model on the finite lattice with periodical boundary 
conditions. It is shown that the partition function of the model on the dual lattice is expressed 
through specific combination of the partition functions on the initial lattice on the torus 
with the corresponding boundary conditions. In section 4 we propose the generalization 
of this relation for the nonhomogeneous case. Here we prove this relation for the weakly 
nonhomogeneous distribution of the coupling constants. In section 5 the exact duality relation 
between two-point order-order and disorder-disorder correlation functions on the torus is 
derived. In section 6 the duality relations for the 2d Ising model with magnetic fields applied 
to the boundaries and the 2d Ising model with free, fixed and mixed boundary conditions are 
obtained. 

2. The model 

The partition function of the 2d Ising model on the torus is 

2 N Z[K}=J2e- ma \ (4) 
M 

-PH[a] = y £*(r)(K 1 (r)V x + K 2 (r)V y )a(r), (5) 

r 

where r = (x, y) denotes the site coordinate on the square lattice of size N = n x m, x = 
1, . . . ,n y = 1, . . . ,m; a(r) = ±1; i^i(r) and K 2 (r) are coupling constants along horisontal 
X and vertical Y axes respectively. The one-step shift operators V x , V y are acting on a(r) 
in the following way 

V x a(r) = a(r + x), V y a(r) = a(r + y), 

V T X = V_„ V T y = V_ y , 

where x, y are unit vectors along horisontal and vertical axes. For periodical boundary 
conditions along X and Y axes 

(vsr = i, (V5r = i, 

and for antiperiodical boundary conditions 

(v£)» = -i, (v;r = -i. 

Here p and a denote periodical and antiperiodical boundary conditions respectively. 

Since we have the four type of the boundary conditions on the torus let us assign the 
corresponding indeces to the partition function: Z^ a '^[K] and, for example, 
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Later on we will consider Z^ a '^[K] as the four-component vector Z[K] with components 
Z b [K], 6 = 1,2,3,4 

Z[K] = (Z ip ' p) , Z M , Z {a ' p) , Z M ). (6) 
We denote site coordinates, functions and functionals on the dual lattice by "tilda" : 

f, a(f), Ki(r), H[a], Z[K], ... 

The site coordinate on the dual lattice coincides with the coordinate of the plaquet center on 
the initial lattice: 

f = r + (x + y)/2. 
Using these denotions, we can write the dual partition function 

2 N Z[K] = Y. e - m \*l 

R 

with the dual Hamiltonian 

-PH[a] =Y,<r(r)(Ki(r)V-x + K 2 (r)V- y )a(r). 

r 

The coupling constans -fQ(r) and Ki(r) are connected by duality condition (3). 

It is known that the partition function of the 2d Ising model can be represented as the sum 
of the functional integrals over the lattice real fermion field [4,5]. For our task it is important 
that this representation can be exactly obtained for the nonhomogeneous model [6-8] : 

Z (P,P) = _ ^_Q(P,P) [X] + QM + Q(a,p) + Q(a,a) ^ ^ 

Here Q^ a '^[K] is the Gaussian functional integral over the four-component Grassmann field : 

qW)[ K ] = (j[ coshK i( r )) I ^exp(5 (a ' /3) [^]), (8) 

r,i 

^=nn#). 

r j=l 

ipj(r) is Grassmann field : 

{^(r),^(r')} = 0. 
In (8) the action has the following form 

= E (£ (0) (^)) + CM Mr))) , (9) 



3 



£<°>ty(r)) = £ MrWiir), 

l<i<j 

£(^)(^(r)) = -ti(r)^(r)V^i(r) + t 2 (r)^ 2 (r)V^ 4 (r), 
tj(r) = tanh.£Q(r). 

Representation (7) was written for periodical boundary conditions. But it is obviously that 
we can write the similar relations for arbitrary combinations (periodical and antiperiodical) 
of the boundary conditions along X and Y axes. For example, Z^ a ^ is distingushed from 
Z(p,p) ^ change of the sign of coupling constants Ki(r) in the last right lattice column : 

Z ( P , P ) _^ z (a, P ) af Ki ( n> y ) _> _ Kl ( n? y ) , 

and therefore 

hin^^-hin.y). (10) 

But transformation (10) is equvalent to the change of the boundary conditions for the one-step 
operator in action (9) 

Therefore, Z^ a ' ,p) is distingushed from Z^ p ' p ^ by the arrangement of signs "±" before terms in 
(7): 

ZM[K\ = ^(q m [K\ + Q M [K] - Q {a ' p) [K] + Q {a ' a) [K]y (11) 

For the rest boundary coditions (Z^ p ' a \ Z^ a ^) we have analogous expressions. Using (6), one 
can write expressions for the partition functions with corresponding boundary conditions on 
the torus in the following compact form: 

Z[K] = RQ[K\. 

where we introduce vector 



and matrix R : 

R 



Q = (Q(P<w,Q(P' a ),Q( a <w,Q( a < a )), 



1 

2 



1 \ 
1 
1 

"I / 



R 2 = 1. 



;i2) 
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Functional integral (8) is expressed through the Pfafhan of the AN x 4iV-dimensional 
antysymmetric matrix D : 



qW)[ K ] = (^J]cosh^(r)jpf(D (a ' /3) ), 



where 



£)(a./3) 





-1 

■l-*i(r)VS 



1 





l + *i(r-x)V^ 1 \ 

1 l + t 2 (r)V% 







V 



-1 -l-t 2 {r-y)Vt y 
and its defines the quadratic form in (9) 



) 



Similar expression one can write for the dual lattice : 



Z[K] = RQ[K], 



where matrix D has the following form 



£)(«./*) 



V 



-1 

-1 -*i(r + x)V c 3 
-1 



1 

-1 

-1 -t 2 (r)Vt 



l + *i(f)V» a 
1 


-1 



1 

1 




(13) 



(14) 



(15) 



ti(r) = tanhi^j(r) 

In conclusion of this section we note that for the finite lattice the Pfafnan of matrix (13) 
or (15) is the finite power polinom of U(r) or ti(r) . It is not hard to calculate Pf(D) in the 
high-temperature limit : 

Pf(£>) = 1, when U(r) = 0, (16) 
and Pi(D) in the low-temperature limit: 



Pf(D) = 1, when U(r) = 0. 



(17) 
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3. The homogeneous case 

In general case for arbitrary sizes m, n and distributions of the coupling constants -fQ(r) 
we can not calculate neither the Pfaffian nor the determinant of matrix D . But in the 
homogeneous case, when 

Ki(r) — Ki — const 

and the matrix D becomes the translation-invariant matrix the determinant of D is calculated 
by means of the Fourier transformation: 

det(D) = 111(1 + tf)(l + t\) - 2ti(l - t\) cosp x - 2t 2 {\ - t\) cospj, 
p 

where U = tanh.£Q and the momentum components p x ,p y (p = (Px,P y )) take the integer and 
half-integer values (in units of 2n/n and 2n/m) for periodical and antiperiodical boundary 
conditions respectively. Using relation 

(Pf(D)) 2 = det(D), (18) 

we obtain from (14), (18) 

Q 2 (K X , K 2 ) = Il( c i c 2 - s i cos ^ - s 2 cospy), (19) 
p 

where 

d = cosh2i^j, Si = sinh2i^"j. 
In similar way one gets for the dual lattice: 

Q 2 (K 1: K 2 ) = n(^ic 2 - h cosp x - s 2 cosp y ). (20) 
p 

Here 

Cj = cosh2i^j, Si = sinh2_ft'j. 
Using (19) and (20), it is not hard to check the following equality 

M'^Q 2 ^, K 2 ) = (hh'r N/2 Q 2 (Ki, K 2 ), (21) 

where Sj and Sj satisfy relations 

Si-S 2 — 1, s 2 -si — 1. 

Note that (21) one can consider as the duality relation for the square of the functional integrals 
which appear in representation (7) for the partition function. Extracting the square root from 
the both parts of (21), we obtain 

{s lS2 y N ' 4 Q^\K u K 2 ) = ±(s 1 s 2 r N ^Q^\K l ,K 2 ). (22) 
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Show that sign "— " in (22) appears only for function Q&' p \ but for the rest components of 
vectors Q and Q we have the sign "+" : 

(aia 2 )- N/4 Q(^i, %) = (8 l82 )- N ^gQ(K u K 2 ), (23) 
where signature matrix g has form 



9 = 



( 


-1 











\ 







1 


















1 







\ 











1 


/ 



Since Pf(D) is the polinom over ti, then from (19), (20) it follows that det(D) is square 
of this polinom. Really in product (19) the multipliers appear by pairs according to the 
momentum components ±p x and rkp y . The exclusion is multipliers corresponding to the 
values p x = p y = it and p x = p y = 0. Let us denote them by q n and q^. 

q n = cic 2 + si + s 2 = (1 + ti + t 2 - t^) 2 cosh 2 K 1 ■ cosh 2 K 2 

Qo — CiC 2 — Si — s 2 = (1 — ti — t 2 — tit 2 ) 2 cosh 2 Ki ■ cosh 2 K 2 = 

= (1 - sxs^/q^. 

Hence it is clear that in functions Q^ 01 '^ (except of Q^^) all multipliers have the constant 
sign in the following range of parameter values : 

si > 0, s 2 > 0, 

and, therefore, these functions do not change the sign. On the other hand at crossing through 
critical line s 1 s 2 = 1 function Q^'^ which contains multiplier (qo) 1 ^ 2 ~ (1 — SiS 2 ) changes the 
sign . Similar results we obtain for dual functions Q^ a '^ . 

Therefore, taking into account the high-temperature and low-temperature limits (16) and 
(17) for the Pfaffian, one gets 

sign(g^(X 1 ,X 2 )) = sign(l - SlS2 ) 

sign(Q^(i? 1 ,i? 2 )) = sign(l - S&) = 
= -sign((Q^\K l ,K 2 )). 

This relation proves (23). 
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Multiplying the right and left parts of (23) on the matrix R (12), we obtain the duality 
relation for the partition functions: 



(s 1 s 2 )~ iV / 4 Z(ir l! ir 2 ) = ( SlS2 )" iV / 4 fZ(^ 1) ir 2 ) ) 



(24) 



f = RgR, f 2 = 1, 
/ 1 1 1 1 \ 
-111-1-1 

T = - . 25 

2 1-11-1 K J 

{ 1 -1 -1 l) 

From (24) it follows that difference between the right and left parts of Kramers- Wannier 
duality relation (1) 

(s 1 s 2 )' N/4 Z^ p \K 1 ,K 2 ) - (s^y^Z^iK^Kz) = 

= (s^r^Q^iKuKi) = -( SlS2 )- JV / 4 g^(X 1 ,X 2 ) (26) 
is equal to zero only on critical line 



sinh 27^ • sinh2iT 2 = 1. 

Moreover outside this line the right and left parts of (26) are compared among themselves at 
increasing the lattice size (m,n — > oo). Therefore duality relation (1) is correct only in the 
thermodynamic limit: 



lim f — In ( 



Z{K U K 2 



lim (UJ Z ^ K - 



m,rwoo \N \ ( Sl S 2 ) N/4 



Since for the finite lattice partition functions Z^ a '^ are analytic functions of Ki (the finite 
power polinom over e ±Ki ) the nonanalytic dependence in the right and left parts of relation 
(24) can appear only from multipliers (s 1 s 2 )~ N ^ 4 ~ ■ Therefore, taking into account the circuit 
rules of branching points corresponding to sinh 1K\ = and sinh 2K 2 = 0, duality relation 
(24) can be continued to the rest ranges of values of the coupling constants: K\ > 0, K 2 < 0; 
Ki < 0, K 2 > 0; K x < 0, K 2 < 0. 
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4. The nonhomogeneous case 



In previous section it was shown that the Kramers- Wannier duality relation (1) (having 
rather symbolic meaning) can be modificated to the exact equality. The KadanofT-Ceva 
anzats (2) also was proved in [3] for special (nonphysical) boundary conditions. Note that the 
thermodynamic limit for the nonhomogeneous Ising model can be defined for specific class of 
the functions determining distribution of the coupling constants -fQ(r). For illustration let us 
consider, for example, the sequence of functions {K^ N \r)} defined in the following way: the 
coupling constants are equal zero on the boundary T of the finite-size clasters: Ki(r e T) = 0. 
In this case at N — > oo we obtain the increasing numbers of the non-interacting finite-size 
clasters. Then this is rather similar on the self-averaging procedure for disodered systems than 
on the usual thermodynamic limit. It means that the duality relation for the nonhomogeneous 
system must be formulated for the finite lattice. 

The covariant notation of exact duality relation (24) for the homogeneous model suggests 
the obvious recipt for the generalization to the nonhomogeneous case. For finite lattice on 
the torus Kadanoff-Ceva anzats (2) must be modified by the following way: 



Unfortunely we can not prove this relation for arbitrary lattice size and distribution of the 
coupling constants. But we checked duality relation (27) for lattices of small size by direct 
calculation on the computer. Moreover the duality relation can be proved for the weakly- 
nonhomogeneous case: 



n(sinh2^(f)) 



1/4 Z[K] = n(sinh2^(r))- 1/4 f Z[K]. 



(27) 



r,i 



Ki(r) = Ki + 5Ki(r), Ki = const, SK^r) < 1, 



when n and m are arbitrary and finite. 
For the first order over SK^r) we have 




-1/4 



r,i 



(28) 



Z[K] = Z(K U K 2 ) 1 + ]T((a(ry(r + x))6K 1 (r) + 



r 



+ {a{r-)a{r- + y))5K 1 {r)) . 



(29) 
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Note that at different K\ and K 2 and (or) n ^ m correlation functions (a(r)a(r + x)) and 
(<j(r)a(r+y)) do not equal among themselves and due to translation invariance do not depend 
from r : 



Z { K 1 ,K 2 )(a { rMr + 1 )) = ^R^§^l 
Z(tf„ K,)(a{r)a(r + y)) = Lr?9^B±. 



(30) 



Taking into account (28)-(30), in the first order over SKi(r) and SKi(r) we obtain for duality 
relation (27): 

-1/4 



[SlS 2 



^ /2[lS-^-Q)E^(r) 



N dK 2 2s 2 



(sis 2 ) 



-1/4 



Z(K U K 2 )+R 



1 <9Q 



Cl 



N dK x 2s! 



Q )E^i(r)+ 



(31) 



,^9^2 2s 2 

For the order zero terms this equality is satisfied according to homogeneous duality relation 
(24). Show that it is satisfied for the linear terms over 5K and 5K terms. From (3) and (23) 
it follows: 



5K x {r) = -^5K 2 (r), 5K 2 (r) = ~5K x {r), 

$2 S\ 







_ d 
-s 2 - 



d 



„ d 

-Si- 



(32) 
(33) 



OK! " z dK 2 dK 2 'dKi 

Substituting (32) and (33) in the left part of (31) and collecting similar terms, one gets 

1 <9Q 



(sis 2 ) 



-N/4 



TZ(K U K 2 ) + Rg 



^Q)£^ 2 (f)+ 



N 8K 2 2s 2 



+ Rg 



1 <9Q 



Cl 



(34) 



MdK! 2si 

Using (25), (from which follows Kg = TR) it is not hard to show that (34) coincides with the 
right part of (31). This proves duality relation (27) in the weakly-nonhomogeneous case for 
arbitrary m and n. 
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In conclusion of this section note that for analysis of the duality relation for correlation 
functions it is convenient to use the other normalization in (27). Using relations 

cosh 2 2ifi(r) _ cosh 2 2K 2 (f) 
sinh2 J ftT 1 (r) ~ s inh2^ 2 (r) ' 

cosh 2 2K 2 (r) _ cosh 2 2% (f) 
sinh2K 2 (r) ~ s inh2X 1 (f) ' 
which follow from (3), and introducing denotions according to [3] 

Y[K] = l[(cosh2K l (r))- 1/2 Z[K}, 

r,i 

Y[K] = l[(cosh2K t (r)r 1/2 Z[K], (35) 



we obtain instead of (27): 



Y[K] = TY[K]. (36) 



5. Duality relation for correlation function 

The duality relation for the nonhomogeneous Ising model is usefull for study the correlation 
function properties. For this aim it is convenient to use the magnetic dislocation representa- 
tion for correlation functions [3] . This representation is based on the obvious equality 

e &-**m ai(T2 = i ai(T2e K™. 

Taking into account that 

(J\O n = (o"lO"2)(o"20"3) • • • (On-lOn), 

one can write 

W M 
where Hamiltonian f3H'[a] of the Ising model with defect differs from f3H[o\ by the change 
of the coupling constants K on K ' = K + in/ 2 on the defect line IV [3] which connect sites 
r and r'\ 



Ki{r) = K, 

K + in/ 2 on the links belonging to path IV 
K on the rest links, 



K[{r) 
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7 is the length of this path ( the number of the "spoilt" bonds). Then, using functionals (35), 
we obtain representation for the two-point correlation function (we omit boundary condition 
indeces) : 

G a (r,r>) = (a(r)a(r')) = Y[K'\/Y[K\. (37) 

In the work [3] it was introduced the correlation function of the disorder parameter //(f). 
This variable characterizes the degree of disorder near the point f on the initial lattice and its 
one can consider as the result of the duality transformation for Ising spin a(r). Correlation 
function (/t(f), fJ>(r')) is determined by means of magnetic dislocation T M : 

G»(r,ri) = (^(?'))=Y[K>yY[K], 

where 

\ —K on the links intersecting of path T M 
| K on the rest linkss. 
The duality relation for correlation functions [3], 

mm) = Hr)a{r')) (38) 

follows from (2) and the transformation of magnetic dislocation T a on the initial lattice to 
magnetic dislocation T M on the dual lattice by means of mapping [3] 

K^r) + m/2 -> K 2 (r) ■ e~™, 
K 2 {r)+nx/2^K 1 {¥)-e" l \ 

which follows from (3). Since duality relation (27) for the finite lattice on the torus differs 
from Kadanoff-Ceva anzats (2) the duality relation for correlation functions on the torus has 
more complicate form. For example, using (36), we obtain for the dual lattice with periodical 
boundary conditions 

Gf p \r,r') = Y M [K"]/Y M [K] = (TY[K'}) M / (TY[K]) M = 

= [Z M G ( P' p) (r,r') + Z (p ' a) G {p ' a) (r,r') + Z^G^' p \r, r') + 

+ Z M Gl a ' a \r,r')]/[Z (p ' p) + Z M + Z {a > p) + Z M }. (39) 

It is not hard to show that (39) pass to (38) only under the following condition: the correlation 
length is smaller of the lattice sizes (this is happen out the scaling domain and at the large m 
and n). Note that duality relation (39) coincides with the relation for correlation functions 
on the torus in the critical point obtained in the paper [9] by means of the quantun conformal 
field theory methods [10]. 
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6. Duality and boundary conditions 

The duality relation for the nonhomogeneous Ising model on the torus allows to obtain 
duality relations for the 2d Ising model with magnetic fields applied to the boundaries and 
for the 2d Ising model with free, fixed and mixed boundary conditions. 

In order to get these relations let us consider the Ising model on the torus with defect 
which is defined by the following distribution of the coupling constants in Hamiltonian (5): 
Ki(r) = K 2 {r) = K on all links of the lattice with the exeption of the following cases - 
Ki{n — l,y) = hi, Ki{n,y) = h 2 , K 2 (n,y) = h, where hi = [3Hi, h 2 = f3H 2 , h = f3H, 
y = 1, ... ,m. Using (3), this defect one can define on the dual lattice by the following way: 
Ki{r) = K 2 (r) = K on all links with the exeption of the following cases - Ki(n — 1, y) — hi, 
Ki{n, y) = h 2 , K 2 (n, y) = h, where hi = (3Hi, h 2 = (3H 2 , h = (3H, y = 1, . . . , m and coupling 
constants h, hi, h 2 and h, hi h 2 are connected by relations: 

sinh2/i • sinh2/i = 1, sinh2/?,i • sinh2/ii = 1, sinh2/i2 • sinh2/i2 = 1. (40) 

Taking limit h — > oo in partition function (2) (for simplicity we consider the ferromagnetic 
model), it is not hard to obtain the partition function of the 2d Ising model with magnetic 
fields Hi and H 2 aplied to the boundaries: 

2Z p {K,h 1 ,h 2 ) = lim (cosh h)~ m Z^ v \K, h x , h 2 ,h) = 

h— >oo 

m m 

2 ^ exp(K ^ a(r)a(r + i) + hi ^ a(n - 1, y) + h 2 ^ a{l,y)), (41) 

[a] r,i y=l y=l 

where in order to sum over spin variables {a(n, y)} we used the following equality: 

m m 

lim (cosh h)" m [] exp(ha(n, y)a(n, y + 1)) = J] S(a(n, y), a(n, y + 1)). (42) 

y=l y=l 

Here on the right-hand side the product of Kronecker's 5-functions is written. From {cr(n, y)} 
this product selects the two spin configurations: all spins are directed up or down. 

Note, that in (41) partition function Z p {hi, h 2 ) is obtained by the limiting procedure with 
corresponding normalization which removes infinite constant. In other cases in consequence 
of the conflict between the product of Kronecker's 5-functions and the boundary conditions 
we can get the zero after taking of the limit, for example, 

lim [2 cosh K)- m Z i - a ' a \K, hi, h 2 ,h) = 0, (43) 

but this does not mean, that Z a (K, hi, h 2 ) = 0. 
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In consequence of (40) for the dual lattice we have h — and 

Z^{K, hi, h 2 ) = \im Z {a ' p \K,h l ,h 2 ,h) = ]T exp(^]T v(r)a(r + i) 

m m 

+ h l J^a{n,y)a{n,y + l) + h 2 J^Z{l,y)d{l,y + l)), (3 = a,p (44) 

y=l y=l 

In result we obtained the partition function of the Ising model on the cylinder with the free 
boundary conditions and the defects on its bases. 

Now, taking limits h — > oo and h — > in (36) and using (41)-(44), it is not hard to get 
the duality relations for the Ising model on the square lattice wrapped on the cylinder with 
magnetic fields applied to its bases 

Z*(K, hi, h 2 ) Z"(K, hi, h 2 ) + ZP(K, hi, -h 2 ) 

(45) 



cosh m 2hi cosh m 2/ i2 (cosh2K) 2 -(™- 3 )] 1/2 [cosh m 2h x cosh m 2/i 2 (cosh 2K)^-V] 1/2 
Z a {K, hi,h 2 ) Z?(K, hi, h 2 ) - ZP(K, hi, -h 2 ) 



• (46) 



cosh m 2hi cosh m 2/i 2 (cosh 2Kf m ^)] 1/2 [cosh m 2hi cosh" 1 2/i 2 (cosh 2Kf m ^)} 1/2 ' 

Here we have lattices with sizes (n — 1) x m and n x m on the right-hand and left-hand sides 
respectively. 

In order to get the partition function of the 2d Ising model on the initial lattice with 
different boundary conditions on the cylinder bases it is necessary to consider in (45)- (46) 
different combinations of limits hi — > oo and hi — > 0, i — 1, 2: 

1) the free boundary conditions - hi — h 2 — 0, which we denote 

Z? 0fi) = 2~ 2m Z a (0,0), a = a,p, 

2) the fixed boundary conditions - hi — > oo: 

Z^ +) = lim(2 cosh hi2 cosh h x y m Z a {h x ,h 2 ), 

Z ( a + _) = lim(2 cosh M cosh h x y m Z a {h x ,-h 2 ), 
3) the mixed boundary conditions - hi — > oo, /i 2 = or /i 2 — > oo, foi = 0: 

Z ( Q +)0) = lim(4cosh/i 1 )- m Z a (/i 1 ,/i 2 ) = lim(4cosh/i 2 )- m Z a (/n,/i 2 ). 

In consequence of (40) the transition to limits hi — > oo and /ij — > 0, i = 1, 2 on the initial 
lattice leads to the following results on the dual lattice: 
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1) h ± ,h 2 ^0, 

\miZ a {KMM) = (2cosh^) 2m Z ( Q 00) , a = a,p (47) 

2) hi — > oo, h 2 — > oo 

lim(2 cosh fci)~ m (2 cosh fc 2 )~ m Z p (fci,^ = 2 [z p (i?, X) + Z P (K, —if)] , (48) 

3) /ii — > oo, /i2 — ► 

lim(2cosh/i 1 )- m Z a (/ il ,/i 2 ) = (2coshK) m 2Z a (K,0), (49) 

Setting in (45), (46) h\ = h 2 = K and respectively h\ = h 2 = K we obtain the first two 
duality relations for the 2d Ising model with the boundary conditions: 

%Q) = Z P (K, K) +JP(K, -K) 
(cosh2iT) m ("- 3 ) (cosh 2K) m ^) 

%o) = Z*(K t K)-Z>(K,-K) 
(cosh 2X) m ("- 3 ) " (cosh 2fsT) m ( n - 3 ) 
Recall that here we have lattices with sizes (n — 1) x m and n x m on the right-hand and 
left-hand sides respectively. 

Now, using (47)- (49) and taking the corresponding limits in (45), (46), it is not hard to 
get the following duality relations: 



(2 cosh Kf™Z\ m _ Zf +|+) + ^ +> _ ) 
(cosh 2iT) m ("- 3 ) (cosh 2K) m ^) ' 

(2 cosh K) 2m %o) _ Zf +>+) ~ Zf + ,-) 



(cosh2iTr( n - 3 ) " (cosh 2K) m ^) ' 
Z p (^, X) + Z p (K, -X) _ Z[ 



7P 

J (o,o) 



(cosh2i^) m ("- 3 ) ~ ~ (cosh 2iT) m ( n - 3 ) ' 
Zp(K,K)-Zp(K,-K) Z? m 



(52) 

(53) 
(54) 
(55) 
(56) 



(cosh2iTr( n - 3 ) " ~ (cosh 2iTT("- 3 ) ' 

Mg^gj) = %+) 
(cosh27Tr( n - 3 ) (cosh2iT) m ( n - 3 )' 

Note that in relations (52), (53) we have lattices with sizes (n — 1) x m and (n — 2) x m 
on the right-hand and left-hand sides respectively. 

Let us make some comments about (53). Using results of exact solution [11] of the 2d 
Ising model with magnetic fields applied to the boundaries, it is not hard to show that 



15 



Z V [K, hi, h 2 ) — Z P (K, hi, —h 2 ) is proportional to sign(hi ■ h 2 ) and therefore at obtaining of 
(46) the contradiction after taking limits hi — > 00, h 2 — > 00 and hi — > 00, /i 2 — > — 00 is not 
appeared. 

In the critical point relations (50)- (55) are reduced to 

z (o,o) - z (+,+) + 

which was obtained in [12] and it follows from [11]. 

7. Conclusion 

From the duality relation for the nonhomogeneous Ising model one can be obtained some 
usefull concequences. Using this relation, one can correctly introduce the "mixed" correlation 
function of type (a /j,a' //) and discover their fermionic content. In principle anzats (27) allows 
to constract the generating functional depending from external currents J(r), J(r) and x( r )> 
where the first two currents are connected with fluctuations of the order and disorder variables 
and the last current generates the fermionic type exitation. 

Undoubtedly the nonhomogeneous duality relation will be usefull for analysis of Ising 
model with random coupling constants. 

However, unfortunaly, we have not proof of this relation for arbitrary distributions of the 
coupling constants and sizes of the lattice. In given paper the duality relation is proved for 
the homogeneous case and in the first order for weakly nonhomogeneous case (it is not hard 
to prove one in the second order). One can prove this relation for case when the some small 
numbers of the coupling constants is chosen arbitrary ones on the background of the rest 
homogeneous coupling constants. The additional argument for correctness of duality relation 
(27) is the direct check of one on the small lattices with sizes m, n — 2, 3. 

We would like thank A. Belavin, VI. Dotsenko, M. Lashkevich for usefull remarks con- 
nected with subject of this paper. 

V.S. thanks Dr. A. Morozov for the hospitality and the exellent conditions at ITEP, where 
this paper has been finished. 
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